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Selected results on Strong and Coulomb-induced correlations from the STAR experiment 
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Using recent high-statistics STAR data from Au+Au and Cu+Cu collisions at full RHIC energy I discuss 
strong and Coulomb-induced final state interaction effects on identical (n — Jl) and non-identical (n — E) particle 
correlations. Analysis of 7t — S correlations reveals the strong and Coulomb-induced FSI effects allowing for 
the first time to estimate space extension of 3t and Z sources and average shift between them. Source imaging 
technique providing clean separation of these effects from effects due to the source function itself is applied to 
one-dimensional relative momentum correlation function of identical pions. For low momentum pions and/or 
non-central collisions large departure from a single-Gaussian shape is observed. 
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I. INTRODUCTION 

Progress in understanding space-time structure of multipar- 
ticle production via femtoscopy technique is currently driven 
by high-statistics data sets accumulated in heavy ion experi- 
ments at RHIC and SPS accelerators Ullilllia]. In partic- 
ular ambitious program of the STAR collaboration at RHIC 
exploiting good particle identification has already provided 
vast variety of femtoscopic measurements in different iden- 
tical and non-identical particle systems some of which were 
actually measured for the first time |4]. This contribution is 
a progress report on two currently pursued STAR femtoscopy 
analyses. Both of them were introduced already at the previ- 
ous WPCF meeting last year in Kromenz |0]- First focuses 
on non-identical particle correlations in rather exotic meson- 
baryon system — |5]. Previous investigations have 
shown that coiTelations among these two charged hadrons re- 
veal not only the Coulomb but also the strong final state inter- 
action (FSI) effects [11, 12]. Order of magnitude difference 
in mass plus AB=1/AS=2 gap in baryon/ strangeness quantum 
numbers makes K — Z system an important tool to study the in- 
terplay between matter flow of on partonic and hadronic level. 
The second analysis exploits particle coiTelations in a more 
conventional system of two identical charged pions (d]. Its 
aim is to understand the geometry of the source. Ultimate and 
rather ambitious aspiration of this project is to extract pion- 
pion scattering lengths. This goal will make heavy ion femto- 
scopic measurements fully competitive to a dedicated particle 
physics experiments trying to extract this important parame- 
ter of the strong interactions |3, 7, 8]. Basic prerequisites for 
such measurements are good knowledge of correlations due to 
quantum statistics and Coulomb interactions. 



II. STUDYING SPACE-TIME STRUCTURE OF 
MULTI-STRANGE BARYON SOURCE VIA 7i S 
CORRELATIONS 

In ultra-relativistic heavy-ion collisions at RHIC hot and 
dense strongly interacting matter is created exhibiting prop- 
erties of de-confined partonic matter jol [Tsll . Almost instan- 
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taneous equilibration of produced matter indicated by recent 
heavy flavor measurements ll9l [l4ll represent one of the great- 
est puzzles coming from the RHIC f^. The early partonic col- 
lectivity also shows up in a subsequent evolution of the system 
leading to strong collective expansion of the bulk matter as 
demonstrated by large values of observed elliptic flow |9, 13]. 
These observations are further substantiated by STAR multi- 
strange baryon measurements showing that E and Q. baryons 
reveal significant amount of elliptic flow which is compara- 
ble to ordinary non-strange baryons flStl . The sizable multi- 
strange baryon elliptic flow which obeys constituent quark 
scaling 1 16] confirms that substantial part of the collective mo- 
tion has indeed developed prior to hadronization. This picture 
is also corroborated by more recent STAR analysis of the cen- 
trality dependence of hyperon yields caiTied out within the 
framework of a thermal model [17]. Observed scaling behav- 
ior of strange baryons is consistent with a scenario of hadron 
formation from constituent quark degrees of freedom through 
quark recombination provided that the coalescence took place 
over a volume that is much larger than the one created in any 
elementary collisions. 

These observations fit nicely into ideal hydro evolution 
starting from the system of de-confined QCD matter. How- 
ever they are also consistent with more realistic hybrid macro- 
scopic/microscopic transport approach |20] which takes into 
account the strength of dissipative effects prevalent in the lat- 
ter hadronic phase of the reaction. The hybrid model calcu- 
lations indicate that at top RHIC energy the hadronic phase 
of the heavy-ion reaction is of significant duration (at least 
10 fm/c) making hadronic freeze-out a continuous process, 
strongly depending on hadron flavor and momenta. In par- 
ticular heavy hadrons, which are quite sensitive to radial flow 
effects, obtain the additional collective push created by res- 
onant (quasi)elastic interactions during that fairly long-lived 
hadronic rescattering stage [21]. 

It is clear that question concerning multi-strange 
baryon decoupling from hot and strongly interacting 
partonic/hadronic system is an interesting one but also not an 
easy one to solve. Could this be provided by the femtoscopy? 
What kind of relevant information can be obtained via low- 
relative-velocity correlations of multi-strange baryons with 
the other hadrons? Since non-identical particle correlations 
are sensitive not only to the extent of the source, but also 
to the average shift in emission time and position among 
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different particle species |19] the answer is affirmative. The 
femtoscopy of non-identical particles was already used to 
show that the average emission points of pions, kaons and 
protons produced in heayy-ion collisions at SPS and RHIC 
are not the same d S [13 I^JH. In hydrodynamically 
inspired blast-wave approach lisll mass ordering of average 
space-time emission points of different particle species natu- 
rally appears due to the transverse expansion of the source. 
This effect increases with a mass difference of the measured 
particle pair Hence studying correlations in the K — Z system 
where the mass difference is really big should provide rather 
sensitive test of the emission asymmetries introduced by the 
transverse expansion of the bulk matter 

Moreover, in addition to the Coulomb interaction studied in 
the K~ K system the small relative momentum — cor- 
relations may provide sufficiently clear signal of the strong 
interaction revealing itself via E*(1530) resonance. Express- 
ing particle momentum in the pair rest frame k* = = — ps 
via pair invariant mass M,t£ and m„ and 
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TABLE I: 200GeV Au+Au, Run IV data set 



Centrality 


No. of S± 


No. of E- 


No. of E+ 


-10% 


1084 X 10^ 


595 X 10^ 


489 X 10^^ 


10 - 40% 


412 X 10^^ 


226 X 10^ 


186 X 10^ 


40 - 80% 


145 X 10^ 


79 X 10^ 


66 X 10^ 



B. Data analysis and corrections 

Event mixing technique was used to obtain uncorrected 
two-particle distribution in pair rest frame. To remove spu- 
rious correlations of non-femtoscopic origin the uncorrelated 
pairs were constructed from events with sufficient proximity 
in primary vertex position along the beam direction, multiplic- 
ity and event plane orientation variables. Pair cuts were used 
to remove effects of track splitting and merging. Resulting 
raw correlation function was then corrected for purity of both 
particle species. The correction was performed individually 
for each bin in k* — (A;*,co.s'9,(p) of the 3-dimensional (3D) 
correlation function as described bellow. 



one expects the E*(1530) peak to show up in the correlation 
function C{k*) at k* 150 MeV/c. Due to its rather long life- 
time Ts*(i53o) = 22/m/c the resonance could be also a rather 
sensitive to the n — Z interaction during long-lived hadronic 
phase. This should be investigated by both the 7t — S fem- 
toscopy as well as via direct measurements of the E*(1530) 
spectra. While first signal of the E* resonance in heavy ion 
colUsion was seen in femtoscopy analysis just two years ago 
||24^ first preliminary STAR measurements concerning the E* 
spectra and their yields were presented only recently at the 
Quark Matter conference in Shanghai this year li25ll . 



A. Data selection 

Though in previous analyses IH [Til the 7t — S 

correlations were studied for two different system dn-Au and 
AuH-Au and also at two different energies in this contribution 
I will concentrate on AuH-Au data at \/Jsnn) = 200 GeV from 
RHIC Run IV only. The data were divided into several cen- 
trality classes. During the run central trigger was used to en- 
hance fraction of 10% most central events. Track-level cuts 
based on dE / dx particle identification in the STAR Time Pro- 
jection Chamber were used. Pion sample momenta p, were 
limited to [0.125, 1.] GeV/c. After the dE /dx cuts the upper 
PT-Mrmt is 0.8GeV/c and 0.6 GeV/c at y = 0. and y = 0.8, re- 
spectively. Charged E were reconstructed topologically in the 
Pt range [0.7, 3.] GeV/c. To increase total number of ana- 
lyzed % — Z pairs in this analysis we have used wider rapidity 
cut then in the previous STAR femtoscopy analyses l26ll . The 
cut \y\ < 0.8 instead of \y\ < 0.5 was employed for both parti- 
cle species. Total number of extracted Z used in this analysis 
are given in Table U 



C. Pair purity analysis 

Pair purity defined as a fraction of primary K — Z pairs 
is calculated as a product of particle purities of both par- 
ticle species. Z purity was obtained from reconstructed Z 
invariant mass plot as a function of transverse momentum. 
Pion purity was estimated using parameter VX of the stan- 
dard parametrization of the identical K — K correlation func- 
tion. The identical pion measurements were performed with 
the same pion cuts as those used in the K~Z analysis. Since 
value of the X parameter is influenced by decays of long-lived 
resonance as well as by non-Gaussian shape of the correlation 
function the pion purity correction can be a significant source 
of the systematic error. 

In order to make contact with previous STAR identical 
pion analyses [26.1 on Fig. 1 we present A;7'-dependence of 
parameters X, Rout, Rside and Riong entering standard out- 
side-long decomposition of 3D correlation function C(q) = 
l+X- exp{-ql,Rl, - q^.,^R^.,^ - ^L/L,,)- Here kr = 
(|pi7'| + |pir|)/2is average transverse momentum of two pi- 
ons. On the same figure the ratio Rom /Rside is plotted too. We 
conclude that the improved cuts used in the present analysis 
do not change the values of extracted parameters but due to 
increased acceptance in rapidity and transverse momenta of 
the pions our analysis covers also region of lower kj- 

We have also studied influence of electrons on purity of 
the pion sample. Exclusion of pions with dE /dx within ±2a 
around the electron band has changed the value of parameter 
X in the kj interval [0.15, 0.25] GeV/c by 50% at maximum. 
However the other parameters characterizing the 3D correla- 
tion function of identical pions (the radii Rout,Rside,Riong) as 
well as the %~Z correlation functions remained unchanged. 
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FIG. 1 : Comparison between parameters of 3D Gaussian fit the 
correlation functions of identical charged pions produced in Au+Au 
collisions at 200 GeV Previous (open stars) and this analysis 
(full triangles). Error bars contain only statistical uncertainties. 



D. Results in ID source size information 

The 7t — E correlation functions C{k*) were analyzed in the 
pair rest frame (k* = Pn = — ps)- The results for unlike sign 
pairs are presented on Fig 2. 
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FIG. 2: The centrality dependence of the correlation function of 
combined sample of unlike-sign It^S^ pairs . 

For all centralities the low k* region is dominated by the 
Coulomb interaction. The strong interaction manifest itself in 
a peak around k* « 150 MeV/c corresponding to S*(1530). 
Peak's centrality dependence clearly shows high sensitivity to 
the source size. Contrary to the Coulomb region the correla- 



tion function in the resonance region does not suffer from the 
low statistics and can thus in principle be used to extract sizes 
of the sources more effectively then in the former case. 



E. Results in 3D asymmetry measurement 

The information about shift in average emission points be- 
tween 71 and E can be extracted from the angular part of the 
3D correlation function C(k*) = C(F,cos9,(p). 
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FIG. 3: Ao()(fc*) and Ai\(k*) coefficients of spherical decomposition 
for combined sample of unlike-sign 71^ pairs from three different 
centrality bins 

The function is binned in k* , cos 9, cp with A^os e = ™d 
= ^ as bin sizes in cos 9 and (|), respectively. After its 
decomposition into spherical harmonics Jz^ 



^ ^ all bins 

A,,„(r) = ^2^ ^ F,,„(9,,(p,)C(r,cos9,-,(p,) (2) 



47t 

symmetry constrains further limit number of relevant compo- 
nents. This is due to the fact that individual coefficients ap- 
pearing in the above decomposition represent different sym- 
metries of the source. Thus for azimuthally symmetric iden- 
tical particle source at midrapidity, only A/,„ with even values 
of I and m do not vanish. On the other hand for non-identical 
particle correlations the coefficients with odd values of I and 
m are allowed. 

For both cases most important coefficient is Aoo(A:*) repre- 
senting angle-averaged correlation function C{k* ). Latter is 



4 



Brazilian Journal of Physics, 2006, in press 



sensitive to the source size only. On the other hand in non- 
identical particle case An (A:*) measures a shift of the average 
emission point in the direction. 

The decomposition of the correlation function into spheri- 
cal harmonics was used to study the centrality dependence of 
asymmetry in emission between pions and E. The results are 
presented on Fig 3. The coefficient A 1 1 which is non-zero in 
all centrality bins clearly confirms that the average space-time 
emission points of pions and E are not the same. 



region dominated by the strong final state interaction the cal- 
culations over predict both, the size and the shift coefficients. 



F. Extracting the source parameters 

For further analysis we have used only the low k* region 
dominated by the Coulomb interaction, excluding thus the re- 
gion around the E*(1530) peak (see Fig. 5). 
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FIG. 4: Comparison o/Aqo (fc* ) and Ai\[k*) coejficients of spherical 
decomposition for combined sample of unlike-sign 71^ pairs from 
10% most central Au+Au collisions with the FSI model predictions. 
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FIG. 5: Comparison o/Aoo(^*) andAii{k*) coefficients of spherical 
decomposition for combined sample of like-sign 71* and unlike- 
sign Tl^E^ pairs from 10% most central Au+Au collisions with the 
FSI model predictions. 



On Fig. 4 the experimental results for the 10% most cen- 
tral, the highest statistics bin are compared to a model calcu- 
lation exploiting strong and Coulomb final state interaction 
lf28ll . Theoretical correlation function was constructed from 
particle emission separation distribution. For this purpose mo- 
menta of real particles were used. The emission coordinates of 
both 71 and E were generated using the blast wave model ifisll . 
This model encompassing correlation between particle mo- 
menta and their space-time coordinates was used with a sin- 
gle set of parameters when generating emission coordinates of 
both K and E. These parameters were obtained from experi- 
mentally measured pion spectra and k — k emission radii. Let 
us note that using the same set of parameters for the E source 
as for the pions implicitly assumes significant transverse flow 
of E. In the Coulomb region theoretical correlation function 
is in qualitative agreement with the data. Moreover, orienta- 
tion of the shift and its magnitude agrees with the scenario in 
which E participates in transverse expansion. However, in the 



The theoretical correlation function was calculated using 
momentum distribution of pairs extracted from the real data. 
The emission coordinates were randomly generated from two- 
parameter source distribution constructed in the following 
way. For both particle species Gaussian shape of their source 
was assumed. The sources were shifted in direction rela- 
tive to each other. This in the pair rest frame can be expressed 
via two parameters R and A„j„ characterizing the pair separa- 
tion distribution Ar*. R representing the width of the Gaussian 
and Aout the shift in the Rout direction. Values of the source 
parameters were then extracted by finding a minimum value 
of between theoretical and real correlation function. Fitting 
was done simultaneously for like and unlike-sign correlation 
functions. 

For most central AuH-Au collision this method yields first 
preliminary values of R— (6.7 ± 1 .0) fm and A„„, — (—5.6 ± 
1.0) fm. The errors are purely statistical. Systematic error 
studies are under way and their values are expected to be of 
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the order of the statistical ones. In our case a negative value of 
the shift means that average emission point of S is positioned 
more to the outside of the whole fire-ball than the average 
emission point of pion as expected in the E flow scenario. 

III. IDENTICAL PION CORRELATIONS 
A. Motivation 

In ultra-relativistic heavy-ion collisions femtoscopy can 
also be employed to extract the information related to the 
strong interaction among the particles llsi fioll . In recent STAR 
analyses the strong FSI was studied via two-ZT^ interferome- 
try I29I] and in proton-lambda correlations fsol] . A model that 
takes the effect of the strong interaction into account has been 
used to fit the measured correlation functions. In the p — A 
and 7? — A correlations, which were actually measured for the 
first time, annihilation channels and/or a negative real part of 
the spin-averaged scattering length was needed in the FSI cal- 
culation to reproduce the measured correlation function. 

At the previous WPCF meeting ambitious proposal to ex- 
ploit the correlations among identical charged pions to ex- 
tract pion-pion scattering lengths was made DSD. Potential for 
such measurement at RHIC and latter on also at the LHC 
stems from the fact that compared to a dedicated particle 
physics experiments measuring scattering lengths parameters 
Qq and like BNL-E865 |31] or Dirac |32] heavy ion ex- 
periments provide much larger number of pion pairs at small 
relative momenta in a single event plus very large data sam- 
ples ( 10^ — 10'^ events). The real challenge when studying 
the strong interaction among identical charged pions then is to 
beat down all systematical errors pertinent to ultra-relativistic 
heavy ion collisions environment. Coulomb interaction, pion 
purity and geometry of the pion source need to be kept under 
control. Varying source size (kr or centrality) may provide a 
good handle on this. 

In particular the bias arising from frequently used Gaussian 
assumption of the source shape needs to be addressed. Using 
high statistics sample of AuH-Au events from STAR experi- 
ment at RHIC highest energy accumulated during the run IV 
it was found that for all kj and centrality bins the Levy stable 
source parametrization does not bring an advantage in describ- 
ing the detail shape of measured three-dimensional correlation 
function [6J. 

Next step in this direction is to exploit model-independent 
imaging technique of Brown and Danielewicz ll33i[34ll and re- 
construct the source itself. This is done via inverting Koonin- 
Pratt equation (see e.g. IjI): 

C{q) -I =4% J Ko{q,r)S{r)r^dr 
Koiq,r) = i| (|<I>(q,r)|2 - l)dicosd^,r) 

1=y-iPi~P2r (3) 

which expresses the ID correlation function C{q) in terms of 
the probability S{r) to emit a pair of particles at a separation 



r in the cm. frame. S{r) is usually called source function. 
All FSI is encoded in the (angle averaged) final state wave 
function <I>(q,r) describing the propagation of the pair from 
a relative separation of r in the pair cm. to the detector with 
relative momentum q. The angle between vectors q and r is 
denoted by 9q r. The procedure to extract S{r) via inversion 
of Eq. 3 is given in Il34ll . 




I I I I I 1 1 1 I I I I I 1 I 

0.02 0.04 0.06 

q._^^ (GeV/c) 



CO 




FIG. 6: ID imaging analysis of isotropic Gaussian source using pion 
momentum spectra from RQMD simulated 10% most central Au+Au 
collisions at y'^Aw =200GeV. Top: original (black filled circles) and 
restored by imaging technique ( red filled circles) correlation function 
C{q). To guide an eye the data points are connected by a smooth 
line of the same color. Middle: ID source function S{r) (red filled 
circles). Full black line represents result of a single-Gaussian fit to 
S(r), dashed line - input Gaussian distribution. Bottom: Atlr^Sij). 



B. Testing the inversion procedure 

To test reproducibility of original source function by the 
imaging procedure we have generated Bose-Einstein (B-E) 
correlated and Coulomb interacting pion pairs. The momen- 
tum spectra of the pions were obtained from 10% most cen- 
tral AuH-Au RQMD generated events at ^/saw =200GeV. The 
pions coordinates were randomly sampled from an isotropic 
3D Gaussian distribution with the radius R=5 fm. The in- 
stantaneous emission 8t=0 of particles was assumed. Using 
pion momenta and their coordinates the B-E correlations and 
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Coulomb FSI were introduced using procedure described in 
lf36ll . Results for correlation and the source functions are dis- 
played on Fig. 6. In order to compare with the input Gaussian 
source the extracted source function was fitted with Gaussian 
distribution again: 



r2 



(4) 



Let us note that the r^-weighting of S{r) appearing in the 
normalization condition of the source function: 



ID imaging analysis of the pion source from Au+Au and 
Cu+Cu collisions 



The same data set of Au+Au events from Run IV used in 
already described 7t — E correlation analysis was also used for 
the reconstruction of the pion source. Second data set em- 
ployed in imaging analysis consists of about 8 million Cu+Cu 
minimum bias events at ^snn =200GeV accumulated dur- 
ing the Run V. Selected results from ID imaging analysis of 
the pion source using centrality selected Au+Au and Cu+Cu 
events are displayed on Fig. 7 and 8, respectively. 



471 / S{r)r^dr - 



(5) 



makes normalization constant X more sensitive to behavior of 
S{r) at r ^ R than one would naively expect. Since latter 
is determined by the imaging from values of the correlation 
function C{q) at very small q fulfilment of the above normal- 
ization condition using data with limited statistics in the low q 
bins may be hard to achieve. This may explain why the values 
of extracted parameters R and X from the single-Gaussian fit 
were found to be smaller than the input ones. While for R the 
discrepancy is about 4% for X it is almost 25%. 
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FIG. 7: ID imaging analysis of identical charged pions from Au+Au 
collisions at y^s^N =200GeV. Results are shown for 10% most cen- 
tral (the left panel) and for peripheral centrality 40-80%) collisions 
(the right panel). Top: measured correlation function C(q) (black 
filled circles ), restored correlation function from imaging technique 
(red filled circles). To guide an eye the data points are connected 
by a smooth line of the same color Middle: ID source function 
S{r) ( red filled circles). Full black line represents result of a double- 
Gaussian fit to S{r). Two dashed lines show contribution of each of 
two Gaussians. Bottom: 4nr^S{r): (red filled circles). Full black 
line represents result of a double-Gaussian fit. 
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FIG. 8: ID imaging analysis of identical charged pions from Cu+Cu 
collisions at y^SiviV =200GeV. Results are shown for 10% most cen- 
tral ( the left panel) and for peripheral ( centrality 40-60%) collisions 
(the right panel). Top: C(q). Middle: ID source function S{r). Bot- 
tom: 47tr^5(r). Labels are same as on Fig.7. 

Compared to model example discussed in the previous 
paragraph the extracted source functions now develop long 
tails 13511 and cannot thus be described by a single Gaussian 
distribution. The only exception are the data from 10% most 
central Au+Au collisions (Fig.7.) To account for observed 
behavior we use simplest extension of Eq.4 and assume that 
the source function obtains contribution from two Gaussians. 
While the first Gaussian g{r,R\) contributing with fraction 
(1 — a) is responsible for the long tails the second one g{r,R2) 
contributing with weight a of width the R2 < Ri describes the 
source function at small r. The double Gaussian distribution: 

Gir,Ri,R2) = {I -a)g{r,Ri) +agir,R2),Ri > Ri (6) 

now seems to describe data quite well. Comparing Au+Au to 
Cu+Cu collisions of the same centrality we conclude that the 
long tails are more pronounced for the smaller system. 

All source functions presented so far were obtained from 
data integrated over the whole range of average pair trans- 
verse momenta kT= [0.125, 0.600] GeV/c). In the following 
I will present results showing how much observed departure 
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from a single Gaussian shape depends on particle transverse 
momenta. 
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FIG. 9: ID imaging analysis of identical charged pions from 10% 
most central Cu+Cu collisions at y^^Aw =200GeV. Results shown 
are for the lowest ( the left panel) and for the highest ( the right panel) 
kj-hin. Top:correlation function C(q). Middle: 1 D source function 
S(r) Bottom: 4ltr^5(r). Labels are same as on Fig.7. 

Such analysis is performed on Fig.9-11. While for Au+Au 
collisions the statistics allowed us to use 9 bins in fcj- for 
Cu+Cu we have used only 4 bins. On Fig. 9 the lowest kT= 
[0.125, 0.250] GeV/c and the highest kT= [0.450, 0.600] 
GeV/c bins from Cu+Cu 10% most central collisions are com- 
pared. We see that the long tail present in the lowest kj bin 
which is characterized by the Gaussion with the width R^l 
fm completely disappears from the highest kj bin. 

Fig. 10 and 11 show A;?- -dependence of the double Gaus- 
sian fit parameters for Cu+Cu and Au+Au collisions at three 
different centralities, respectively. We observe that with in- 
creasing kr contribution of the second Gaussian is less and 
less prominent. This behavior is more pronounced in Cu+Cu 
collisions. Let us note that long tails observed in the low-A:^ 
bins, then may indicate an important role played by the pi- 
ons from long-lived resonance decays. An interesting obser- 
vation which may have some relevance to observed behavior 
regarding a random walk nature of the particle freeze-out in 
coordinate space was made by T. Csorgo at this meeting iIstIi . 
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FIG. 10: kj -dependence of the source parameters of identical 
charged pions from Cu+Cu collisions at i/JjvA' =200 GeV taken at 
three different centralities: 0-10% - full red circles, 10-40% - full 
blue squares, 40-60% - full black triangles. 
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FIG. II: kj -dependence of the source parameters of identical 
charged pions from Au+Au collisions at ^/smn =200 GeV taken at 
three different centralities: 0-10% - full red circles, 10-40% - full 
blue squares, 40-80% - full black triangles. 



High-Statistics STAR data from Au+Au and Cu+Cu colli- 
sions taken at full RHIC energy and three different centrali- 
ties were used to discuss recent progress in identical (n — K) 
and non-identical (n — S) particle femtoscopy. In 7t — E sys- 
tem the strong and Coulomb-induced FSI effects are observed 
making it possible for the first time to estimate the average 



shift and width between % and E source. ID imaging of iden- 
tical pion source reveals significant departure from a single 
Gaussian shape. Observed long tails which could be fairly 
well described by allowing for the second Gaussian are more 
pronounced for the sources produced at large impact parame- 
ters. The effect is stronger for Cu+Cu then for Au+Au col- 
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lisions. For all centralities and both colliding systems the 
highest departure from single Gaussian shape is observed for 
sources emitting particle with a small average pair transverse 
momenta kj- For highest analyzed kj the source is to a good 
approximation a single Gaussian. 
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